Surface modes and photonic modes in Casimir 
calculations for a compact cylinder 



V V Nesterenko 

Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 
Dubna 141 980, Russia 

E-mail: nestr@theor.jinr.ru 

Abstract. A rigorous formulation of the problem of calculating the electromagnetic 
vacuum energy of an infinite dielectric cylinder is discussed. It is shown that the 
physically relevant spectrum of electromagnetic excitations includes the surface modes 
and photonic modes. The mathematical procedure of summing over this spectrum is 
proposed, and the transition to imaginary frequencies is accomplished. As a result, it 
is justified the imaginary-frequency representation for the vacuum energy which has 
been used in previous Casimir studies for this configuration. 
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1. Introduction 

The notion of the elementary excitation spectrum is of paramount importance in 
all condensed matter physics [H [2]. The excitations of different type result, as 
a rule, in different physical consequences. Therefore it is of a certain interest 
for theoretical and experimental investigations of the Casimir effect to answer the 
question: the electromagnetic oscillations of what type are considered in the problem at 
hand [21 HJ El El E E] • Having elucidated this point one can hope to link, in a transparent 
way, the Casimir force with actual physical properties of the material boundaries. 
However it is not easy to answer this question even when the Casimir force is calculated 
by making use of the familiar Lifshitz formula [9l[T0l[lT] . A rather complicated derivation 
of this formula [12] in the original papers initiated its obtaining anew by making use 
mainly of the mode summation method p^l [T3l [T^| [15^ [T6l [T71 [18^ [T9] . 

The Casimir calculations for nonflat boundaries turned out to be much more 
involved in comparison with those for planes [20]. Especially complicated calculations 
have been done for a circular cylinder |22l 1231 12H I2S1 12SI I2ZI I2SI HSl IBOl I3H 
EH [351 EE! E3 EE]- It was also unclear electromagnetic excitations of what kind have 
been taken into account in these studies [39] . 

The present work seeks to present a consistent derivation of the formula for 
the vacuum energy of electromagnetic field connected with a material cylinder by 
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summing explicitly the contributions to this energy given by different branches of the 
electromagnetic spectrum in this problem. 

The layout of the paper is as follows. In Sec. 2 the spectral problem generated 
by the Maxwell equations for a compact infinite cylinder is formulated rigorously and 
the physically relevant spectrum of electromagnetic excitations for this configuration is 
determined. It is shown that this spectrum includes surface modes (bound states) and 
photonic modes. In Sec. 3 the summation over this spectrum is accomplished by making 
use of the spectral density when accounting of the photonic (continues) branch of the 
spectrum. In Conclusion (Sec. 4) the meaning of obtained results are discussed briefly. 



2. Physical spectrum of electromagnetic excitations for a cylinder 

In the source-free case the general solution to Maxwell equations can be represented in 
terms of two independent Hertz vectors [4"0] 



e = v x v x n' + i/i^ v x n", (i) 

H = - % e - V x n' + V x V x n" . (2) 

Here II' is the electric Hertz vector, II" is the magnetic Hertz vector, c is the velocity of 
light in vacuum, and the Gauss units are used. The Hertz vectors obey the Helmholtz 
vector equation 

(v 2 + P)n = o, (3) 

where the wave number k is given by 

J 1 

k 2 =£[! — . (4) 

cr 

On the other hand, it is known that the general solution to Maxwell without sources 
equations can be derived from two scalar functions which may be chosen in deferent 
ways J4TJH2]. F° r the configuration with cylindrical symmetry the role of these functions 
can play the axial components of the electric (II') and magnetic (II") Hertz vectors. 
The rest components of II' and II" are zero in this case. As a result, the Helmholtz 
vector equation ([3]) reduces to the scalar Helmholtz equations for 11^, = n' and II" = n" 

LU 2 \ 

a + eii— J n = o, n = n', n" (5) 

with the following general solutions 

n'= £ a n f™(r)e ihz+in9 , (6) 

n=0,±l,±2,... 

n"= Yl b n f™(r)e lhz+me . (7) 

n=0,±l,±2,... 

The cylindrical coordinates (r, 9, z) are used and the z axis coincides with the axis of a 
circular infinite cylinder of radius a. The medium inside the cylinder has the permittivity 
E\ and permeability These quantities outside the cylinder acquire the values 62 and 
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fj,2, respectively. We assume for definiteness that E\\l\ > £2^2- The wave vector along 
the z axis is denoted by h. The amplitudes a n and b n for the solutions inside the cylinder 
will be denoted by a l n and b l n , respectively, and in the same way for solutions outside 
the cylinder we introduce the amplitudes a e n and b e n . 

The functions and f™{ r ) m the general solutions (jHJ) and (J7J) obey the 

radial wave equation 



dr 2 r dr \ r 2 / 

Inside the cylinder we put 

fn(r) = J n (Xir), n = 0,l,..., 0<r<a, (9) 

where J n {z) is the Bessel function, Ai = \Jk\ — h 2 , k\ = u 2 / c\ = e\^\UJ 2 jc 2 . Outside 
the cylinder we consider first 'outgoing' waves 

/„(r) = #«(A 2 r), n = 0,l,..., r>a, (10) 



where H^\z) is the Hankel function of the first kind, A2 = y/k^ — h 2 , k\ = oo 2 jc\ = 
e 2 ^2 uj 2 /c 2 . 

In the radial solutions (j2J) and (ITDl the sign of X 2 = k 2 — h 2 , s — 1, 2 is not fixed 
yet. Thus, in our consideration the solutions 

/nO) = /„(A» for r <a (11) 

and 

Ur) = K n {\ 2 r) for r>a (12) 

are also admissible. Here X 2 = h 2 — k 2 , s = 1,2, ^(z) = i~ n J n {iz) and K n {z) = 
i n+1 ^H n (iz) are the modified Bessel functions [13]. 

On the cylinder surface the matching conditions should be satisfied. These 
conditions require the continuity of tangential components of fields E and H when 
crossing cylinder surface 

discont (E||) = 0, discont (Hy) = 0. (13) 

The matching conditions give rise to the frequency equation determining admissible 
values of the spectral parameter uj in the boundary value problem under consideration: 

Cc> 2 a 4 / . J' , H' \ ( J' H' 



2 y 7 H I \ J H 

1 \ ,J n 1J -n J \ °n ±± n 

n 2 h 2 a 2 



£iX 2 -r — £2X1-^- ] ( ^1X2-^ — ^Xi-pp- 

Is 

0, // 0.1.2 (11) 



xfx 2 . 

In this equation 



•2 "1 2 

—(exfii - e 2 ^2) 
c z 



J n = J n (X ia ), H n = H^(X 2 a), X s = + v / lu 2 /c 2 s - h 2 , s = l,2, (15) 

the prime on the functions J n and H n means differentiation with respect to their 
arguments, C\ and c 2 are the velocities of light inside and outside the cylinder, 
respectively, c 2 > c\. 
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The roots of (I14p are important in radio-engineering when developing the 
radio dielectric waveguides [HI H3 06j 071 SHI HH1 [50] and in fiber optics (optical 
waveguides [511 E2]). The results of investigation of the frequency equation ( fl4l) 
determining the spectrum in the problem at hand can be summarized in the following 
way. All the real roots of this equation lie in the interval 

c\h < to < c 2 h . (16) 

These roots make up two discrete sequences. In the interval (fTBl) the frequency equation 
ffl~4l) can be rewritten in the form 



wV ( - J' n x K'\ ( - J' K' n 



n 2 h 2 a 2 



\ 2 \ 



2 



uj 2 



— (ei/ii - E2H2) 



C; 



2 



0, n = 0,1,2,... (17) 



with the notation K n = K n (X 2 a). 

When frequency uj equals the real roots of equation (Tl7|) . located in the interval 
( ITBl . the 'outgoing' waves (fTUl) become the functions decaying in the radial direction 
f fT2|) (surface or evanescent waves). These eigenmodes describe the propagation of 
electromagnetic waves along the cylinder (waveguide solutions). The radial functions 
( TTTj) are not realized in the problem under consideration. 

Now we address the complex roots of the frequency equation (fT4l) . In the strip of 
the complex frequency plane 

< Re uj < c 2 h (18) 

there are no complex roots of fTl4|) with Im u ^ 0. 
In the semi-plane 

Reu; > c 2 h (19) 

for sure there are complex roots of fTT4|) with Im w ^ 0. Indeed, in the domain (fT9j) the 
left-hand side of (1141) is a complex function of the complex variable uj. The complex 
eigenfrequencies of a dielectric cylinder lead to leaky (radiating) modes. It is clear that 
these modes cannot carry the electromagnetic energy along the cylinder. For us it is 
important that the modes with complex uj (quasi-normal modes [39]) do not satisfy 
standard completeness condition and as a result they cannot be used for quantization 
of electromagnetic field in the problem at hand. 

In order to get rid of the complex eigenfrequencies and consequently to escape leaky 
or radiating modes we shall consider, outside the cylinder, the scattering states instead 
of outgoing waves. The scattering solutions to Maxwell equations can be derived from 
outgoing solutions by making use of the substitutions 

<#«(A 2 r) -> a+H+(X 2 r) + a-H~(X 2 r) 

blH£\\ 2 r)^b + n H+{\ 2 r)+b-H-{\ 2 r)- 

For simplicity in (120]) the notations 

H+ = Hj», H~ = H^ (21) 
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are introduced. 

As a result, for a given n and h we have 6 amplitudes a l n ,b l n ,a^,a~,b^,b~. The 
matching conditions at the cylinder surface lead to 4 linear homogenous equations for 
these amplitudes. Hence no restrictions arise here for the spectral parameter u 2 /c 2 . 

Eliminating in these equations the amplitudes a l n and b l n we are left with 2 equations 
in 4 amplitudes and b^ 



where 



= " h ( I A l 



a \ Af 



r±' 



The S matrix in this problem 



= UAs/ff - e, A, JffJ ) . 1 23 i 




6+ , ■ , ,- , (24) 

obeys obviously the following matrix equation 

K-S = K + , (25) 

and 

Het K + 

det 5 = P^- . (26) 
det K~ 

By a direct calculation one can easily show that det K~ coincides (up to unimportant 
multiplier) with the left-hand side of the frequency equation (fl4"|) . Thus this equation 
can be rewritten in the form 

det K~ — 0. (27) 



Surprisingly formulas ([26]) and (1271) for the S matrix were not known in the literature 
devoted to the electromagnetic scattering by a cylinder. 

Summarizing we infer that the spectrum of electromagnetic oscillations in the 
problem under study consists of discrete values u> na : c\h < u na < cih corresponding 
to the surface modes and a continuous branch of the spectrum with real positive u: 
c 2 h < uj < oo. In mathematical scattering theory [53J it is proved that the bound states 
and scattering states form together a complete set. 



Surface and photonic modes in Casimir calculations 



6 



3. Summation over the spectrum and transition to imaginary frequencies 

Now we address the calculation of the vacuum energy in the problem at hand proceeding 
from the standard mode-by-mode summation 



E c = - ^2u q = — ^ J^Ungjh) + / uAp n (u,h)duj 

1 {q} J -°° n=0 r» Joih, 



(28) 



where the prime over the sum sign means that the term with n = is taken with 
the weight |. The first term in square brackets is responsible for the surface waves 
contribution and the second one describes the contribution of the photonic modes. The 
latter contribution is represented by making use of the respective spectral shift function 
Ap [39 J. The rigorous mathematical scattering theory gives the following expression for 
the spectral density shift 

Ap(k) = p(k) -po(k) = ^-^ trlnS(fc) = ^^lndetS(fc). (29) 

Here p(k) is the density of states for a given potential (or for a given boundary conditions 
in the case of compound media) and Po(k) is the spectral density in the respective free 
spectral problem (for vanishing potential or for homogeneous unbounded space). It is 
obvious that in the Casimir calculations one has to use just Ap(k) subtracting at this 
point the so-called Minkowski space-time contribution to the vacuum energy. 

In the case of scalar scattering problem the Jost functions f(k) and f(—k), the 
scattering matrix S(k), and the phase shift S(k) are related by the formula 

S(k) = e 2 *' = . (30) 

Substitution of (130]) into (128"]) gives more familiar formula for spectral density [51] 

^^fH^- (31) 

In the problem under consideration the TE and TM modes do not decouple. 
Therefore we are dealing here with the matrix (2 x 2) scattering problem and we must 
use the spectral density defined by (129]) . 

The contribution of the surface modes in ( 128]) can be represented by the counter 
intergral 



(32) 



where F n (u) is the left hand side of (]T7j) . This equation was written for real u. However, 
in the counter integral (132]) an analytical continuation of this function to the complex 
frequency plane should be used. It can be done immediately in terms of det K + (lower 
semi-plane us) and deti^~ (upper semi-plane u). After that we can use for both terms 
in (|28]) the counter integral representations with the counters C + and C_, respectively. 
The counter C_ starts at ioo and goes along the positive imaginary axis to the origin 
and after that it goes along the positive real semi-axis to infinity. The counter C + is 
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obtained by refection of C_ to the upper semi-plane u. As a result we arrive at the 
following imaginary frequency representation of the vacuum energy in the problem at 
hand 



where F n (uj,h) is the left hand side of the frequency equation (|14p . It is this 
representation that has been used in the Casimir calculations for a material cylinder. 

4. Conclusion 

We have shown that in the case of a material cylinder there are two types of 
electromagnetic excitations which are physically relevant: i) surface modes and ii) 
photonic modes. A consistent transition to imaginary frequencies requires the both 
branches of the spectrum are to be taken into account. The contribution to the 
Casimir energy due to the surface modes and photonic modes can be separated only in 
terms of real frequencies. Upon transition to imaginary frequencies these contributions 
are indivisible. Presented consideration justifies rigorously the imaginary frequency 
representation for the Casimir energy of a compact infinite cylinder that has been used 
in many previous papers dealing with investigation of this energy. 

It is worth noting that the mathematical consideration presented here is completely 
applicable to the Lifshitz configuration, namely, to an infinite dielectric plate placed in 
vacuum (dielectric films). 
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